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' Abstract 

I The continuous single-facility min-sum Weber location problem based upon the lift 

' metric is investigated. An effective algorithm is developed for its solution. Implementation 

for both the discrete and continuous location problems is developed in the programming 
package Mathematica. 

AMS Subj. Class.: 90B80, 90B85 
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^ ■ 1 Introduction 

' Location problems represent a special class of optimization tasks, where coordinate of locations 

and distances between them are main parameters. In the general case, the task of location 
problem is to define positions of some new facilities from the actual space in which are already 
placed some other relevant objects (points). New facilities are centers that provide services 
and called suppliers. Existing facilities are the service users or clients, and called customers. 
Location problems occur frequently in real life. Many systems in the public and private sectors 
are characterized by facilities that provide homogeneous services at their locations to a given set 
of fixed points or customers. Examples of such facilities include warehouse location, positioning 
a computer and communication units, locating hospitals, police stations, locating fire stations 
in a city, locating base stations in wireless networks. 

Different classifications of the location problems are known. The classification scheme from 
|14j assumes five positions. 

In the present article we pay attention to the selection of the distance function as the 
characterization criterion of the location problem. The distance between two points is the 
length of the shortest path connecting them. The metric by which the (generalized) distance 
between two points is measured may be different in various instances [3]. In the calculating 
of distance between two points, the most common distance metrics in a continuous space are 
those known as the class of Ip distance metrics, primarily rectangular (/i), Euclidean {I2) and 
Chebyshev (Zoo) metric. Detailed explanation of various metrics one can find in Dictionary 
of distances [6]. Many factors affect on the process of metrics choosing. The most important 
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factor is the nature of the problem. For example, if it is possible to move rectilinearly between 
two points, the distance between them is exactly given by the Euclidean (or straight-line 
distance) metric. On the other hand, in the cities where streets intersect under the right angle 
mainly, the distance between two points will be the best approximated using the rectangular 
metric (also known as the Manhattan, "city block" distance, the right-angle distance metric or 
taxicab distance). Measures of distances in chess are a characteristic example. The distance 
between squares on the chessboard for rooks is measured in Manhattan distance; kings and 
queens use Chebyshev distance, and bishops use the Manhattan distance. 

We emphasize the next main contribution of our paper. 

The Weber location problem (also called the Fermat-Weber problem) is a basic model in 
the location theory which has received significant attention in the scientific literature. For 
a detailed review see, for example, [32]. The paper [23] investigated a reformulation of the 
unconstrained form of the classical Weber problem into an unconstrained minimum norm 
problem. The classical Weber problem is established with the Euclidean norm underlying in 
the definition of the distance function. But, other measures, principally Ip norms, also play 
an important role in the theory and practice of location problems. The norms are arbitrary, 
in general. The most popular method to solve the Weber problem problem with Euclidean 
distances is given by a one-point iterative procedure which was first proposed by Weiszfeld |31j . 
The procedure is readily generalized to Ip distances (see, for example, [19], Ch. 2). Solution of 
the continuous Weber problem in li distance is described in [7]. The three-dimensional Fermat- 
Weber facility location problem with Tchebychev distance is investigated in [25]. The Weber 
location problem with squared Euclidean distances is considered in [7j; the same problem 
under the assumption that the weights are selected from a given set of intervals at any point, 
is studied in [9]. 

The Ip norms have received the most attention from location analysts. But, many other 
types of distances have been exploited in the facility location problem. A review of exploited 
metrics is presented in [7]: 

- central metrics [26], 

- distance functions based on altered norms [201 [2T] , 

- weighted one- infinity norms [30j . 

- mixed norms |15j . 

- block and round norms [27] . 

- mixed gauges [TO] . 

- asymptotic distances [16] , 

- weighted sums of order p [21 [29] . 

In the present article we solve the Weber problem in the plane, under the assumption that 
the distance is measured by the lift metric. 

The paper is organized as follows. Some basic definitions and algorithms are restated in 
the second section. In the third section we present an effective algorithm for the solution of 
the single-facility continuous Weber problem, assuming that distances are measured by the lift 
metric. 
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2 Preliminaries 



The lift metric or the raspberry picker metric in the plane 72.^ is defined by 



L{A,B) 



\x 



+ \x 



X2 



^2 



(2.1) 



where A{xf, x^) and B{xf jxf) are given points. It can be defined as the minimum Euchdean 
length of all admissible connecting curves between two given points, where a curve is called 
admissible if it consists of only segments of straight lines parallel to rr-axis, and of segments of 
y-axis (see, for example [6]). Therefore, under the assumption X2 7^ the distance between 
two points A and B in the lift metric equals the sum of lengths AA' , A'B' and B'B, where 
A' and B' are orthogonal projections of the points A and B to the y-axis, respectively (Figure 
l,Left). In the opposite the distance between A and B is simply the length of 

the segment AB (Figure 1, Right). 




A'= B' 



Figure 1. Left) The case ^ Right) The case = x^ 

This distance is appropriate for usage in cities which have one main street (corresponding 
to the y-axis), and the other side streets are normal to it (Figure 2). We are observed that in 
the main city of Zakynthos island in Greek-Zakynthos, the streets are deployed on this way. 
Similar situation also occurs in tier buildings where the lift (in the role of y-axis) connects 
tiers. 

The 2-dimensional continuous Weber location problem can be briefly restated as follows 
(see, for example [SIE2])- Let m demand centers ^1, . . . , A^ be given in the plane 7^^ (locations 
of given customers), where ^4(05^,02), i = l,...,m. It is necessarily to find a new point 
X{xi,X2) € "R? which has minimal sum of weighted distances with respect to given points. 
Therefore, one needs to solve the unconstrained optimization problem (single-facility min-sum 
problem), where it is necessary to minimize the sum 



min f{X) 



i=l 



Wi ■ d{Ai,X). 



(2.2) 



The real quantity Wi is a positive weighted coefficient of the point Ai. Essentially, the weight 
Wi converts the distance d{Ai, B^) into a cost of serving the demand of customer i considerate 
to A;th offered facility location. 

For the sake of completeness, we restate well-known method which gives solution of the 
Weber problem (|2.2p in the case when the underlying distance function is defined by the h 
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metric, which means 

li{Ai,X) = \xi -a\\ + \x2 - 41- 
Then the goal function f{X) divides in two sums: 

m m 

fii^i) = ^'Wi\xi - a\\, f2{x2) = ^Wi\x2 - al\. 

i=l 1=1 

Therefore, the initial optimization problem splits into two independent optimization tasks. It 
is sufficient to consider the optimization problem 

m 

min f{x) = y Wi\x — a*|. (2-3) 

i=l 

We restate well-known algorithm for solving the problem (j2.3p (one can find it in \23\ 128]). 
The first step in this procedure is optional, but it can be used to accelerate the remaining two 
steps. 

Algorithm 2.1 Solve the single-facility continuous Weber location problem. 

Require: Real quantities a^, . . . , a™. 
1: Step /. For each subset of identical elements a*^ = a*^ = • • • = a*-* perform the following 

activities: put Wi-^ = wi^ + Wi^ • • • + wi^ , eliminate multiple elements a*^ , . . . , a*J as well as 

corresponding weights Wi^ , . . . , Wi. and later perform appropriate shifting of the indices of 

residual elements and their weights Wj. 
2: Step //. If Step / is applied, denote by q the cardinal number of different elements in the 

set {a}, . . . , a™}; otherwise, use q = m. Sort coordinates o*, i = 1, . . . , g in non-descending 

order. Let the sorted sequence of coordinates is 

a^' <a^' <■■■ < a'''. 

Rearrange the weighting coefficients {wi, . . . , Wq} {wy, . . . , } applying identical re- 
placements on the weights. For the sake of simplicity, let us denote partial sums of the 
array {wr,...,Wq>} by 

k 

S2[0]=0, S2[k]=Y,w^', l<k<q. (2.4) 
1=1 

3: Step ///. There are two possibilities, denoted by Pi and P2- 
Pi . If the condition 

S2[k* - 1] < ^S2[q] < S2[k*] (2.5) 

is satisfied for some k* € {!,..., g}, then we end the algorithm without any possible 
solution. Indeed, formal solution x = is eliminated according to assumption ()3.9p in 
this case. 

P2 . If the condition 

\s2[q]=S2m (2.6) 

holds for some k* G {1, . . . then the solution is multiple, i.e. the searched coordinate 
X2 can to have any value from the interval {a^ ,a'^*^^). As agreed, we use the midpoint 
value X = (a^* + a^*~^^)/2. 
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3 Continuous Weber problem and lift metric 

In the sequel we solve the single-facility min-sum Weber problem ()2.2p applying the lift metric 
(j2.ip . Therefore, the distance function is defined by 

d{A,x) = LiA,x)f^ M+lx^-a^l + Kl, Z7t ' ^ = i'---'-- M 

Two major steps (denoted as Step 1 and Step 2) are separated in our algorithm, as in 
the following. 

Step 1. Generate the list X of permissible solutions of the problem. Its initial value is the 
empty set X = 0. Two different procedures are separated during the construction of the set X, 
in accordance with the definition (13. ip . 

Procedure 1. Let us consider the quotient set S/= of the set S = {ag, . . . , a™} into the 
equivalence classes S = {Si = [02], . . . ,Sd = [a-f]}^ where the class Sj contains elements from 
S whose values are 02 ■ For each j = 1, . . . ,d seek the second coordinate of the optimal point 
X{xi,X2) in the form 

X2 € Sj 4^ X2 = ■ (3-2) 

So, as we know value of the coordinate X2 {x2 = 0,2), it is necessary to determine value of 
the coordinate xi. Denote by Qj the set of indices corresponding to points whose second 
coordinates are contained in the set Sj. According to (13. ip and (13. 2p . the objective function 
f{X) consists of two separated sums 

m m 

f{X) = (|xi| + \x2 - 4| + |ai|) + Wi\xi-al-i\ (3.3) 

Taking into account X2 = and grouping the first term in the first sum with the second sum, 
we obtain 

m m 

f{X) =Y,Wi {\a^ -all + 14 1) +Y^wi\xi - (3.4) 

where 

0, m^i- ^ ^ ^ 

As the first sum in the expression (|3.4p is constant, the problem is reduced on determining the 
minimum of the function 

m 
i=l 

We apply Algorithm 12.11 in adapted form for our specific situation (j3.5p . (|3.6p . That process 
consists of three major steps. 

Step /. For each subset of identical elements a^^*^^ = 4^*^^ = . . . = af^^^^ perform the 

following: put Wi-^ = Wi^ + Wi^ + • • • + Wi^ , eliminate multiple elements well 
as their weights , . . . , Wi- and then perform appropriate renumeration of the indices of the 

remainder elements af^''^ and their weights wj. 
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Step //. If Step / is applied, denote by p the cardinal number of different elements in the set 
{a^^^\ . . . , a^^"^^}; otherwise, use p = m. Sort coordinates a^*-^\ . . . , a^^^^ in non-descending 
array. Furthermore we suppose that the ordered sequence is 

Ol < Oi < • • • < ^1 • 

Rearrange the corresponding weighting coefficients wi^ . . . ,Wp analogously in the sequence 

Wli,W2', . . . ,Wp>. 

For the sake of simplicity, let us denote partial sums of the array {wy , . . . , Wpi } by 



5i[0]=0, 5i[A:] = ^Wj/, 1 < A; <p. 



1=1 

Step ///. There are two possible cases capable to produce permissible minimizers for /i, 
denoted by Ci and C2. 
Ci . If the inequalities 

S,[k'-l]<^Si[p]<Si[k'], (3.7) 

are satisfied for some k' € {1, . . . ,p}, then the searched coordinate is xi = a^' . Later, we use 
X{xi,X2) as the possible optimal point: X = XL) {{ai ,02)}. 
C2 ■ If the condition 

]^SM = Si[k'] (3.8) 

is satisfied for some A;' G {1, . . . then the solution is multiple, i.e. the searched coordinate 
xi can to have any value from the interval [a^ , all ~^^]. In our implementation we use the value 
xi = [a^ + all '^^)/2. Thus, we found the additional possible solution of the starting problem 
([22]), which implies: X = X U {((of + of +^)/2, 4^)}. 
Procedure 2. Compute X2 under the assumption 

X2 ^ S ^ X2 for each i € {1, . . . , m} (3-9) 

(under the assumptions opposite with respect to (|3.2p ). the function /(^) is reduced to 

m 

f{X) = Y,w^ + \X2 - 41 + |al|) . (3.10) 
i=l 

It is necessary to minimize that function. Since the third term ^i)j|a|| in the function f{X) 
defined in ()3.10p is constant, one needs to minimize the next two objectives: 



mm fi{xi) = y^Wi\xi\ (3.11) 

i=l 
m 

mm f2{x2) = ^Wi\x2 - aW- (3.12) 



1 = 1 



Thus, solving the problem (|3.10p with two variables was reduced to solving two independent 
tasks of unconstrained optimization (j3.1ip and (j3.12p with one variable (xi and X2, respec- 
tively) . 
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Solution of the optimization problem (j3.1ip is evidently xi = 0. The optimization problem 
(I3.12P is the classical continuous Weber location model with underlying li metric, where only 
the additional assumption (j3.9p is imposed. Therefore, in order to find optimal value for X2 it 
suffices to apply Algorithm 12.11 assuming that the input sequence is = a\, . . . , a™ = and 
taking into account conditions ()3.9p . 

Step I. For each subset of identical elements = = • • • = perform the following 
activities: put Wi-^ = Wi^ + Wi^ + • • • + Wi- , eliminate multiple elements a]^ ^ . . . , as well 
as corresponding weights Wi^ , • • • , Wi- and later perform appropriate shifting of the indices of 
residual elements a!^ and their weights Wj. 

Step II. If Step / is applied, denote by q the cardinal number of different elements in the 
set {a2, . . . ,a2*}i otherwise, use q = m. Sort coordinates a^, i = I, ■ ■ ■ ,q in non-descending 
order. Let the sorted sequence of coordinates is 

^' 2' a' 
0"l < 02 5; ■ ■ ■ < O2 • 

Rearrange the weighting coefficients . . . , Wq) \wv ^ ■ ■ ■ , Wq'} applying identical replace- 
ments on the weights. Subsequently, generate the partial sums 5'2[i],i = 0,. . . ,q of the array 
{wi',. . . , Wqi } as in (j2.4p . 

Step III. There are two possibilities, denoted by Pi and P2. 

Pi. If the condition (|2.5p is satisfied for some k* € {l,...,q}, then the algorithm is 
finished without any solution. Indeed, the formal solution X2 = eliminated according to 

assumption (13. 9p . actual for this case. 

P2. If the condition (12. 6p holds for some k* € {1, . . . ,q}, then the solution is multiple, i.e. 
the searched coordinate X2 can to have any value from the interval (a2*,02 ~^^)- the 
midpoint value X2 = (02 +'^2 ~^^)/'^^ that the possible optimal point is X(0, X2). Place the 
point X at the end of the list Xhy X = XU {(0, (af + af +^)/2)}. 

Step 2. Thus, we got one or more permissible solutions of the starting problem ()2.2p . For 
all obtained values X from X we determine the values of the function f{X) defined in (I3.10p . 
Solution of the Weber problem will be the point X* (x'^, X2) for which the function f{X) has a 
minimal value. Actually in this step we are solving generated discrete location problem, where 
the set X contains in advance defined feasible locations of the supplier. 

Let Xi, . . . , Xr be r locations on which it is possible to set a new desired object (supplier). 
The sum of weighted distances from the permissible location X^, A; € {1, . . . , r} of the supplier 
to the customers is equal to 

m 

Wk = Y,Wi-L{Ai,Xk). (3.13) 

i=l 

The task is to determine the location B^* for which the sum of weighted distances is minimal, 
i.e. 

Wfc. = min {Wk\ l<k<r}. 
In accordance with the previous considerations, we state the following general algorithm. 

Algorithm 3.2 Solution of the single-facility min-sum Weber problem in the lift metric. 

Require: List Ip = {(a}, af), . . . , (a^, a^)} and the list of weights It = {wi, . . . , Wm}- 
1: Step 1: Form the quotient set of 5 = {02, . . . , o™} in the form {Si, . . . , S^}, where each 
equivalence class Sj contains identical elements from S with the value • 
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2: Step 1: Generate the list X applying the procedure included into the possibilities Ci,C2 
(included in Procedure 1.) to all distinctive values of the set S, i.e. using X2 = a-^i 
j = l,...,d. 

3: Step 1: Extend the list X applying the method defined in the case P2 (included in Pro- 
cedure 2.). 

4: Step 2: Solve the discrete location problem using given locations Ip, discrete set X of 
possible solutions and the weights It. 

Example 3.1. Solve Weber problem using the specified algorithm with the next data: 

^1(4,4), lui =4, ^2(3,1), W2 = 1, ^3(6,4), 1^3 = 2, ^4(6,2), it;4 = 3. 

We have S = {4, 1, 4, 2}. The quotient set of S is defined as Si = [4], 5*2 ~ [1], S3 = [2]. Therefore, it is necessary 
to consider three possibilities for the cases Ci and C2. 

1. Let be X2 — 4. Then the function f{X) has the following form: 

f{x) = W2(|2;i| + \x2 - all + l^il) + TO4(|a;i| + \x2 - at\ + \at\) + wi\xi - ai\ + wa\xi - a?|. 

According to the constant value a;2 = 4 of the coordinate X2, the function / just depend on xi, so we can 
consider the next function 

fi{xi) = wi|a::i - a\\ + W2\xi — 0| + W3\xi ~ al\ + W4\xi - 0| 
4 

w,\xi - aj^ ^1, 

i = l 

where af = = 4, af = a? = 6, = af ^ = 0. 

Let us sort the coordinates af'''^ — >■ a\ and rearrange corresponding weights Wi — > , using the same 
replacements: 



Table 1. 



coordinates (a\ ) 








4 


6 


weights (wii) 


1 


3 


4 


2 


k 


1 


2 


3 


4 




1 


4 


8 


10 



We firstly assume that Step / is omitted. According to ^ X]i=i ™i' = 5, the condition 

SiW -1] < ^Si[4] < Si[k'] 

is satisfied for k' — 3, so xi — af = 4. Therefore, one possible solution is Xi(4, 4). 
In the case when Step I is applied, data from Table 1 reduce to 



Table 2. 



coordinates {a\ ) 





4 


6 


weights (wii) 


4 


4 


2 


k 


1 


2 


3 


Silk] = 


4 


8 


10 



Then conditions p.7|) are satisfied for A: = 2, so that the same possible solution is generated. 

The list of permissible solutions is now equal to X = {^i}- 

2. In this case it is assumed X2 = 1. Now the function fi{xi) looks like: 

fi(xi) = wi\xi - 0| + ii;2!a:i - 3| + W3\xi - 0| + TO4|a;i - 0|. 
On the similar procedure as in the case 1. we get the table: 
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Table 3. 



coordinates (a\ ) 











3 


weights (wii) 


4 


2 


3 


1 


k 


1 


2 


3 


4 


Silk] = ELi™'' 


4 


6 


9 


10 



Inequalities (|3.7|) are valid for k' = 2, so Xi = of = 0, i.e. we got the second possible solution ^2(0, 1). 
In the case when Step I is applied, Table 3 transforms to 



Table 4. 



coordinates {a{ ) 





3 


weights (wii) 


9 


1 


k 


1 


2 


= ELi'"«' 


9 


10 



Condition (|3.7|) holds for fc = 1, so that X2 is again the second eventual solution. 
We have X = {Xi,X2}. 

3. Let us now start from the assumption X2 = 2. 

/i(xi) = wi\xi — 0| + W2\xi — 0| + walxi - 0] + w;4|a::i — 2|. 

The relational table is: 



Table 5. 



coordinates (a\ ) 











2 


weights (wii) 


4 


1 


2 


3 


k 


1 


2 


3 


4 


Silk] = ELi ™»' 


4 


5 


7 


10 



Since the equality of the form p.Sp are satisfied for k' — 2, the solution xi is from the interval [af , af ] = [0, 0], 
i.e. xi = 0, which implies X3{0, 2). 

The list X is expanded: 3£ = {^1, X2, X3}. 

Let us observe that Step I transforms Table 5 into the next table 



Table 6. 



coordinates {a\ ) 





2 


weights (wii) 


5 


3 


k 


1 


2 


Silk] = ELi "^i' 


7 


10 



Now, condition (|3.7p is satisfied for fc = 1, so that ^3(0, 2) is possible optimal point. 

4. Let be X2 7^ 1, 2, 4. In this case we take xi — and then seek for the minimum of the function 

4 

/2(x-2) = ^Wi\x2 - all 

i=l 

Let sort the coordinates 02 — >■ 02 and perform analogous rearrangement Wi Wii . The relational table is 



Table 7. 



coordinates [a\ ) 


1 


2 


4 


4 


weights {Wii) 


1 


3 


4 


2 


k 


1 


2 


3 


4 


S2ik] = ELi™»' 


1 


4 


8 


10 
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Inequalities of the form ((23} hold for k* = 3. We stop algorithm. This case has no solution, since the assumption 
X2 7^ 4 is made. 

Let us mention that Step I gives the next Table 8. 

Table 8. 



coordinates (02 ) 


1 


2 


4 


weights (wii) 


1 


3 


6 


k 


1 


2 


3 


S2[k] = T,Li'^i' 


1 


4 


10 



Therefore, the conclusion is the same as from Table 7. 

At the end, in order to solve Step 2 of Algorithm 13.21 we must compute and compare the values of the 
function / at each point Xi, i — 1,2, 3. We get 

/(Xi) = 50, /(X2) = 55, /(X3) = 62. 

Therefore, the solution of the Weber problem will be the point in which the function / has a minimal value, i.e. 
X* =Xi = Ai = (4,4). 

4 Conclusion and future work 

Our paper is the first attempt to solve the discrete and the single-facihty min-sum continuous 
location problem with the lift metric as the measure of distances. 

A couple of variants and extensions of continuous location problems have been investi- 
gated in literature. Let us mention main between them. More complex problems include the 
placement of multiple facilities. Problems with barriers are the subject in [5l 113 ^ \T7 \ I18j. The 
location of undesirable (obnoxious) facilities requires to maximize minimum distances (see, 
e.g., [U ttH [El [22] . Location models with both desirable and undesirable facilities have been 
analyzed in [1] . It seems interesting to investigate these extensions in the sense of the lift met- 
ric or in the more general nonconvex case, where the shortest length of arc is used as distance 
instead of a particular metrics. 

References 

[1] J. Brimberg, A. Mehrez, Multi-facility location using a maximin criterion and rectangular distances, 
Location Science 2 (1994), 11-19. 

[2] J. Brimberg, R.F. Love, Properties of ordinary and weighted sums of order p used for distance 
estimation. Recherche Operationnelle, 29 (1995), 59-72. 

[3] R. Chen, Location problems with costs being sums of powers of Euclidean distances, Comput. & 
Ops. Res. 11 (1984), 285-294. 

[4] P.C. Chen, P. Hansen, B. Jauniard, H. Tuy, Weber's problem with attraction and repulsion. Journal 
of Regional Science 32 (1992), 467-486. 

[5] P.M. Bearing, K. Klamroth, R. Seargs JR., Planar location problems with block distance and Bar- 
riers, Ann. Oper. Res. 136 (2005), 117-143. 

[6] E. Deza and M.M. Deza, Dictionary of Distances, Elsevier, Boston, 2006. 

[7] Z. Drezner and H. Hawacher, Facility location: applications and theory, Springer- Verlag, Berlin, 
Heidelberg, 2004. 

[8] Z. Drezner, A note on the Weber location problem, Ann. Oper. Res. 40 (1992), 153-161. 

[9] Z. Drezner, C.H. Scott, On the feasible set for the squared Euclidean Weber problem and applica- 
tions, Europ. J. Oper. Res. 118 (1999), 620-630. 
[10] R. Durier, C. Michelot, Geometrical properties of the Fermat- Weber problem, J. Math. Anal. Appl. 
117 (1985), 506-528. 



Single-facility Weher Location Problem based on the Lift Metric 11 



[11] E. Erkut, S. Neuman, Analytical models for locating undesirable facilities, Europ. J. Oper.Res. 40 
(1989), 275-291. 

[12] F. Follcrt, E. Schomcr, J. Scllcn, Subquadratic Algorithms for the Weighted Maximum Facility 

Location Problem, http : / / www . nipi-sb . mpg . dc / ^ schocmer/publicat ions / WM, 1995. 
[13] H.W. Hamacher, S. Nickel, Combinatorial algorithms for some 1-facility median problems in the 

plane, Europ. J. Oper. Res. 79 (1994), 340-351. 
[14] H.W. Hamacher, S. Nickel, Classification of location models. Location Science 6 (1998), 229-242. 
[15] P. Hansen, J. Perreur, J.F. Thisse, Location theory, dominance and convexity: Some further results, 

Oper. Res. 28 (1980), 1285-1295. 
[16] M.J. Hodgson, R.T. Wong, J. Honsakcr, The P-Centroid Problem on an Inclined Plane, Operations 

Research 35 (1987), 221-233. 
[17] B. Kafer, S. Nickel, 2001. Error bounds for the approximate solution of restricted planar location 

problems, Europ. J. Oper. Res. 135 (2001), 67-85. 
[18] K. Klamroth, Planar Weber location problems with line barriers. Optimization 49 (2001), 517-527. 
[19] R.F. Love, J.G. Morris, G.O. Wesolowsky, Facilities Location: Models and Methods, North-Holland, 

New York, 1988. 

[20] R.F. Love, J.G. Morris, Computation procedure for the exact solution of location- allocation problems 
with rectangular distances. Naval Research Logistics Quarterly 22, (1975), 441-453. 

[21] R.F. Love, W.G. Truscott, J. Walker, Terminal Location Problem: A Case Study Supporting the 
Status Quo, J. Opl Res. Soc. 36 (1985), 131-136. 

[22] E. Melachrinoudis, An efficient computational procedure for the rectilinear maximin location prob- 
lem. Transportation Science 22 (1998), 217-223. 

[23] N. Mladenovic, Continual location problems, Matematicki institut, SANU, Beograd, 2004, In 
Serbian. 

[24] LA. Osinuga and O.N. Bamigbola, On the Minimum Norm Solution to Weber problem, SANSA 
Conference Proceedings (2007), 27-30. 

[25] M. Parthasarthy, T. Hale, J. Blackhurst and M. Frank, The three-dimensional Fermat-Weber prob- 
lem with Tchebychev distances, Advances Modeling and Optimization (2006), 65-71. 

[26] J. Perreur, J. Thisse, Central Metrics and Optimal Location, Journal of Regional Science 14 (1974), 
411-421. 

[27] J.F. Thisse, J.E. Ward, R.E. Wendell, Some Properties of Location Problems with Block and Round 
Norms, Operations Research, 32 (1984), 1309-1327. 

[28] V.A. Trubin, Effective algorithm for the Weber problem with a rectangular metric. Cybernetics 
and Systems Analysis, 14(6), DOI:10.1007/BF01070282, Translated from Kibernetika, No. 6, pp. 
67-70, November-December, 1978. 

[29] H. Uster, Weighted sum of order p and minisum location models, Ph.D. Thesis, McMaster Univer- 
sity, Canada, 1999. 

[30] J. Ward, R. Wendell, A New Norm for Measuring Distance Which Yields Linear Location Problems, 

Operations Research, 28 (1980), 836-844. 
[31] E. Weiszfeld, Sur le point pour lequel la somme des distances de n points donnes est minimum, 

Tohoku Math. J. 43 (1937), 355-386. 
[32] G. Wesolowsky, The Weber problem: History and perspectives, Location Science 1 (1993), 5-23. 



